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ABSTRACT 

Wc present U BVIc CCD photometry of the young open cluster Stoek 8 with 
the aim of studying its basic properties such as the amount of interstehar 
extinction, distance, age, stellar contents and initial mass function (IMF). We 
also studied the star formation scenario in this region. From optical data, the 
radius of the cluster is found to be ^ 6' (~ 3.6 pc) and the reddening within the 
cluster region varies from E{B — y) = 0.40 to 0.60 mag. The cluster is located 
at a distance of 2.05 ± 0.10 kpc. Using Ha slitless spectroscopy and 2MASS 
NIR data we identified Ha emission and NIR excess young stellar objects 
(YSOs), respectively. From their locations in the colour-magnitude diagrams, 
majority of them seem to have ages between 1 to 5 Myr. The spread in their 
ages indicate a possible non-coeval star formation in the cluster. Massive stars 
in the cluster region reveal an average age of ^ 2 Myr. In the cluster region 
(r ^ 6') the slope of the mass function (MF), F, in the mass range ^ 1.0 ^ 
M /Mq < 13.4 can be represented by a power law having a slope of —1.38 ± 
0.12, which agrees well with Salpeter value (-1.35). In the mass range 0.3 ^ 
M/Mq < 1.0, the MF is also found to follow a power law with a shallower 
slope of r = — 0.58±0.23 indicating a break in the slope of the IMF at ^ IMq. 
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The slope of the K-hand himinosity function for the cluster (r ^ 6') is found 
to be 0.31 ± 0.02, which is smaller than the average value (~ 0.4) obtained 
for embedded star clusters. 

A significant number of YSOs are distributed along a Nebulous Stream 
towards the east side of the cluster. A small cluster is embedded in the Neb- 
ulous Stream. The YSOs lying in the Nebulous Stream and in the embedded 
cluster are found to be younger than the stars in the cluster Stock 8. The radio 
continuum, MSX, IRAS mid- and far-infrared maps and the ratio of [S II] /Ha 
intensities indicate that the eastern region of Stock 8 is ionization bounded 
whereas the western region is density bounded. The morphology seems to in- 
dicate that the ionization/ shock front caused by the ionizing sources located 
in the Stock 8 region and westwards of Stock 8 has not reached the Nebulous 
Stream. It appears that star formation activity in the Nebulous Stream and 
embedded cluster may be independent from that of Stock 8. 

Key words: stars: formation — stars: luminosity function, mass function — 
stars: pre— main— sequence — open clusters and associations: individual: Stock 
8. 

1 INTRODUCTION 

The study of star formation process and stellar evolution constitute one of the basic problems 
in astrophysics. Since most of the stars tend to form in clusters or groups, star clusters 
are useful objects to study the star formation process. The initial mass function (IMF), 
defined as the distribution of stellar masses at the time of birth is one of the basic tools for 
understanding the formation and evolution of stellar systems. Star forming regions (SFRs) 
and young star clusters have proven to be the ideal laboratories for studying the form of 
IMF and its variation through space and time. Since the young clusters (age < 10 Myr) 
are assumed to be less affected by dynamical evolution, their mass function (MF) can be 
considered as the IMF. However, a recent study by Kroupa (2007) argues that even in the 
youngest clusters, it is difficult to trace the IMF, as clusters evolve rapidly and therefore 
eject a fraction of their members even at a very young age. 

There are many unsolved issues concerning the universality of the stellar IMF. The 
theoretical expectation is that the IMF of a cluster should depend on the location, size, 
metallicity, density of the star forming environment and other conditions such as temperature 
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or pressure (Zinnecker 1986; Larson 1992; Price & Podsiadlowski 1995). But, convincing 
evidence for variation in the stellar IMF has not yet been found observationally (Scalo 1998; 
Kroupa 2001, 2002). Massey et al. (1995) concluded that the slope of the IMF at the higher 
mass range is universal, with a value of r= —1.1 ± 0.1 irrespective of the variations in 
metallicity. Luhman et al. (2000) also found no systematic variations in IMF of the nearby 
SFRs (such as p Oph, IC 348, Trapezium) although there is a large difference in stellar 
densities. On the other hand Hillenbrand (1997) found a few SFRs with unusual IMFs. 

The lack of strong evidence for IMF variations suggests that its fundamental form may 
be universal. Consequently, local conditions may not play any significant role in the star 
formation process. However, the influence of star formation history of young clusters on the 
form of IMF is still an open issue. High-mass stars have strong influence on their surroundings 
and can significantly affect the formation of low mass stars. Recently, a relatively large 
number of low mass stars have been detected in a few OB associations, e.g. Upper Scorpio, 
the (J and A Ori regions (Preibisch & Zinnecker 1999, Dolan & Mathieu 2002). Since this 
realization, surveys have demonstrated that the IMFs must essentially be the same in all 
star forming regions. The apparent difference is due mainly to the inherent low percentage 
of high mass stars and the incomplete survey of low mass stars in high-mass star forming 
regions (e.g. Preibisch & Zinneker 1999, Hillenbrand 1997, Massey et al. 1995). Recently 
Parker & Goodwin (2007) have studied the origin of 0-stars and their effect on low-mass 
star formation. Advancement in detectors along with various surveys such as the 2MASS, 
DENIS and available archived data from ISO and Spitzer space telescopes have permitted 
detailed studies of low-mass stellar population in regions of high mass star formation. 

With the aim of understanding the star formation process and IMF in and around young 
star cluster regions containing OB stars, we selected young cluster Stock 8 located within 
the HII region of IC 417 (Sh2-234) in the Auriga constellation of the Perseus arm. No CCD 
photometric observations have been carried out for this region so far. Mayer & Macak (1971) 
have carried out photoelectric photometry of 11 bright stars in the region and reported a 
distance of 2.96 kpc for the cluster. Using the photographic study of stars brighter than V 
= 16 mag, Malysheva (1990) found cluster age as ~ 12 Myr with an angular diameter of 20' 
and a distance of 1.9 kpc. 

The present study is an attempt to understand the structure, star formation history, 
pre-main-sequence (PMS) population and the form of IMF of Stock 8 using the UBVIc 
CCD photometry data. In Sections 2 and 3, we describe the observations, data reductions 
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and archived data used in the present work. Sections 4 to 8 describe various resuhs. The 
star formation scenario in the IC 417 (Sh2-234) region is described in Sections 9 to 12. The 
results are summarized in Section 13. 



2 OBSERVATIONS AND DATA REDUCTIONS 

The 30 X 30 arcmin^ (~ 17.9 x 17.9 pc^) area containing the central region of Stock 8, 
reproduced from the DSS2-R band image, is shown in Fig. [H As evident from the figure, 
Stock 8 is embedded in an HII region. To the west of it, there are several OB stars which 
appear to form an OB association. To the east of Stock 8, we find a very peculiar chord-like 
nebulosity which we term as Nebulous Stream. A small cluster is also embedded within the 
Nebulous Stream. The small embedded cluster and the Nebulous Stream are discussed in 
Sections 9 and 11, respectively. We describe below the observations made to carry out a 
detailed study of the region. 

2.1 Optical CCD Photometry 

The CCD U BVIc observations were carried out using the 2048 x 2048 pixeP CCD camera 
on the 105-cm Schmidt telescope of the KISO Observatory, Japan on November 4, 2004. 
At the Schmidt focus (f/3.1) each pixel of 25 fim corresponds to 1.5 arcsec and the entire 
chip covers a field of ~ 50 x 50 arcmin^ on the sky. The read-out noise and gain of the 
CCD are 23.2 e~ and 3.4 e~/ADU, respectively. The FWHM of the star images was ~ 3 
arcsec. A number of bias and dome flat frames were also taken during the observing runs. 
We took short and long exposures in all the filters to avoid saturation of bright stars. SA 92 
field of Landolt (1992) was observed to standardize the field on the same night. The log of 
observations is tabulated in Table [H 

The CCD frames were bias-subtracted and flat-field corrected in the standard manner 



using various tasks available under IRAl]^. Aperture photometry was done for the standard 
stars of SA 92 field to estimate the atmospheric extinction and to calibrate the field. Fol- 
lowing calibration equations are derived using a least-square linear regression: 



{U -B) = (0.887 ± 0.025)(m -h)- (2.973 ± 0.010) 
[B-V) = (1.233 ± 0.009)(6 - v) + (0.482 ± 0.010) 

^ IRAF is distributed by National Optical Astronomy Observatories, USA 
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Table 1. Log of observations 
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^ 105-cm Schmidt Telescope, Kiso, Japan 
^ 104-cm Sampurnanand Telescope, ARIES, Naini Tal 
2-m Himalayan Chandra Telescope, lAO, Hanle 



(V -I) = (0.847 ± 0.008)(t; - i) + (0.758 ± 0.006) 
V = v- (0.095 ± 0m8){V -I)- (2.578 ± 0.007) 



where u, b, v, i are the instrument magnitudes corrected for the atmospheric extinctions 
and U,B,V,I are the standard magnitudes. Fig. [2] shows the standardization residuals, A, 
between standard and transformed V magnitudes, {U — B), {B — V) and (V" — /) colours of 
standard stars as a function of V mag. The standard deviations in AV, A(f/ — 5), A(i? — V) 
and A(V" - /) are 0.021, 0.045, 0.023, 0.015 mag, respectively. 

Different frames of cluster region of same exposure time and filters were averaged. Pho- 
tometry of cleaned frames was carried out using the DAOPHOT-II (Stetson 1987) profile 
fitting software. The magnitudes of bright stars were taken from short exposure frames 
whereas that of faint stars were taken from deep exposures. Profile-fitting photometry gives 
the error in magnitude determination, the goodness of the fit parameter, x, which is a mea- 
sure of the average rms deviation to the PSF fit normalized to the expected errors. It also 
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Figure 1. The 30 X 30 arcmin'^ DSS2-R band image of the region around Stock 8. The area marked with thin lines is the 
region covered by shtless spectroscopic observations and the area marked with thick lines is the region for which deep images 
are taken in V and Ic bands. The circle represents the location of the small embedded cluster (see Section 9). Locations of the 
Nebulous Stream (see Section 11) and Stock 8 are also marked. The abscissa and the ordinates are for J2000 epoch. 



gives a shape parameter, Sharpness, which measures how well the PSF fits the object. The 
photometric errors in magnitude, colours and image parameters as a function of V magnitude 
are shown in Fig. [31 These parameters were used to reject poor measurements. 



Deep images of two central regions of the cluster and a nearby field (02000 = 05^29™'42'^; 
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Figure 2. Residuals between standard and transformed magnitude and colours of standard stars plotted against the Landolt 
standard V magnitude. Theses observations are from Kiso Schmidt telescope. The error bars represent combined errors of 
Landolt (1992) and present measurements. 

(52000 = +34°42'13") were also obtained in V and Ic filters on September 29, 2006 using 
a 1024 X 1024 pixeP CCD camera mounted at the f/13 Cassegrain focus of the 104-cm 
Sampurnanand Telescope (ST) of the Aryabhatta Research Institute of Observational sci- 
encES (ARIES), Naini Tal, India. Log of the observations is given in Table [1] In this set 
up, each pixel of the CCD corresponds to 0.37 arcsec and the entire chip covers a field of 
~ 13 X 13 arcmin^ on the sky. To improve the signal to noise ratio, the observations were 
carried out in binning mode of 2 x 2 pixels. The FWHM of the star images was ~ 1.8 arcsec. 
The observed region is marked by thick lines in Fig. [TJ The secondary standards from the 
KISO observations were used to calibrate the data taken at ARIES. The typical DAOPHOT 
photometric errors at brighter end {V ~ 15) are of the order of ~ 0.01 mag, whereas the 
errors increase towards the fainter end (~ 0.05 at ^ ~ 21). The combined photometric data 
of KISO Schmidt telescope and ST in an area of 40' x 40' along with positions of the stars 
are given in a table, which is only available in electronic form as part of the online material 



from http://www.hlackwell-synergy.com Format of the table is shown in Table [2l 
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Figure 3. The photometric errors in V, (U — B), {B — V), {V — I), image parameters x 'ind Sharpness as a fmiction of V 
magnitude for Kiso observations. 



Table 2. UBVIc photometric data of the sample stars. The complete table is available in electronic form only. 



star a(2ooo) '5{2000) V {U - B) (B-V) (V - I) 

no (h:m:s) (d:m:s) 



1 


05:27:48.31 


+34:21:26.6 


11.210 


0.032 


0.746 


0.724 


2 


05:28:01.22 


+34:27:00.7 


11.498 


0.129 


0.174 


0.281 


3 


05:27:43.61 


+34:21:24.8 


11.271 


-0.481 


0.417 


0.617 


4 


05:28:50.26 


+34:19:23.2 


11.507 


0.015 


0.174 


0.298 


5 


05:28:04.46 


+34:29:20.8 


11.710 


-0.474 


0.298 


0.494 



To study the luminosity function (LF)/ MF it is necessary to take into account the 
incompleteness in the observed data that may occur for various reasons (e.g. crowding of the 
stars). A quantitative evaluation of the completeness of the photometric data with respect 
to the brightness and the position on a given frame is necessary to convert the observed 
luminosity function (LF) to a true LF. We used the ADDSTAR routine of DAOPHOT II to 
determine the completeness factor (CF). The procedure has been outlined in detail in our 
earlier work (Pandey et al. 2001). Briefly, we randomly added artificial stars to both V and 
/ images in such a way that they have similar geometrical locations but differ in I brightness 
according to mean {V — I) colour (~ L5) of the data sample. The luminosity distribution 
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Table 3. Completeness Factor (CF) of photometric data in the cluster and field regions. 



V range 


Stock8 




Field region 


l^mag; 


r ^ 3' 


3' < r sC 6' 




10 - 11 


1.00 


1.00 


1.00 


11 - 12 


1.00 


1.00 


1.00 


13 - 14 


1.00 


1.00 


1.00 


14 - 15 


1.00 


1.00 


1.00 


15 - 16 


1.00 


1.00 


1.00 


16 - 17 


0.98 


0.97 


0.98 


17 - 18 


0.95 


0.95 


0.97 


18 - 19 


0.93 


0.93 


0.96 


19 - 20 


0.88 


0.90 


0.92 


20 - 21 


0.88 


0.86 


0.86 


21 - 22 


0.74 


0.75 


0.76 



of artificial stars is chosen in such a way that more stars are inserted towards the fainter 
magnitude bins. The frames are reduced using the same procedure used for the original 
frame. The ratio of the number of stars recovered to those added in each magnitude interval 
gives the CF as a function of magnitude. The minimum value of the CF of the pair (i.e. V 
and / band observations) for the two sub-regions, given in Table [3|, is used to correct the data 
for incompleteness. The incompleteness of the data increases with increasing magnitude as 
expected, however it does not depend on the area significantly. 

2.2 Slitless Spectroscopy, Ha and [S II] imaging 

The cluster was observed in slitless mode with grism as the dispersing element using the 
Himalaya Faint Object Spectrograph Camera (HFOSC) on the 2-m Himalayan Chandra 
Telescope (HCT) of Indian Astronomical Observatory (lAO), Hanle, India. The central 
2K X 2K pixels of 2K x AK CCD were used for data acquisition. The pixel size is 15 /im 
with an image scale of 0".296/pixel and it covers an area of 10x10 arcmin^ on the sky. 
A combination of Ha broad-band filter (Ha-Br; 6100 - 6740 A) and Grism 5 were used 
without any slit. The resolution of Grism 5 is 870. The regions where slitless spectroscopic 
observations were carried out are identified with thin lines in Fig. [H Emission-line stars with 
enhancement over the continuum at Ha wavelength are visually identified. Multiple frames 
were taken to confirm the presence of Ha emitting sources. Positions of the Ha emission 
stars are given in Table H Ha-Br and [S II] (A = 6724 A, AA ~ 100 A) images of two central 
regions of the cluster were also acquired on September 24, 2006 using the HFOSC at the 
2-m HCT. Log of the observations is given in Table [H 



10 Jessy Jose et al. 

Table 4. Coordinates and photometric data of detected Ha emission stars 
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3 ARCHIVAL DATA 

Near-infrared (NIR) JHKg data for point sources within a radius of 20' around Stock 8 have 
been obtained from the Two Micron All Sky Survey (2MASS) Point Source Catalogue (PSC) 
(Cutri et al. 2003). To improve the photometric accuracy, we used photometric quality flag 
(ph_qual = AAA) which gives a S/N ^ 10 and a photometric uncertainty < 0.10 mag. This 
selection criterian ensures best-quality detection in terms of photometry and astrometry (cf. 
Lee et al. 2005). The JHKg data were transformed from 2MASS system to CIT system 
using the relations given on the 2MASS web site. 

The Midcourse Space experiment {MSX) surveyed the Galactic plane in four mid-infrared 
bands - A (8.28 /im), C (12.13 /im), D (14.65 /xm) and E (21.34 /xm) at a spatial resolution 
of ~ 18" (Price et al. 2001). Two of these bands (A and C) with \{AX) corresponding to 8.28 
/im (3.36 /im) and 12.13 /im (1.71 /im) include several polycyclic aromatic hydrocarbons 
(PAHs) features at 6.2, 7.7, 8.7, 11.3, and 12.7 /im. MSX images in these four bands around 
the cluster region were used to study the emission from the PAHs. 

The data from the IRAS survey around the cluster region in the four bands (12, 25, 
60, 100 /im) were HIRES processed (Aumann et al. 1990) to obtain high angular resolution 

^ http:/ /www. astro. caltech.edu/~jmc/2mass/v3/transformations/ 
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maps. These maps were used to study the spatial distribution of dust colour temperature 
and optical depth. 

4 RADIAL STELLAR SURFACE DENSITY AND CLUSTER SIZE 

The radius of a cluster is one of the important parameters used to study the dynamical 
state of the cluster. We used star count technique to determine the statistical properties 
of clusters with respect to the surrounding stellar background, to study the surface density 
distribution of stars in the cluster region and to derive the radius of the cluster. 

To determine the cluster center, we used the stellar density distribution of stars having 
V ^ 18 mag in a ± 75 pixel wide strip along both X and Y directions around an eye 
estimated center. The point of maximum density obtained by fitting a Gaussian curve is 
considered as the center of the cluster. The coordinates of the cluster center are found to 
be aaooo = 05'^28'"07^5 ± P.O; (52ooo = +34°25'42" ± 15". The 2MASS NIR data yields the 
cluster center ~ 7" away southwards of the optical center. However, this difference is within 
the uncertainty and therefore can be considered identical. Henceforth, we have used the 
optical center as the center of the cluster. 

To investigate the radial structure of the cluster, we derived the radial density profile 
(RDP) using the Kiso Schmidt observations for stars brighter than ^ = 18 mag by dividing 
star counts inside the concentric annuli of 1' wide around the cluster center by the respective 
annulus area. The densities thus obtained are plotted as a function of radius in Fig. HI where 
one arcmin at the distance of cluster (2.05 kpc, cf. Sec 6) corresponds to ~ 0.6 pc. The 
upper and the lower panels show the RDPs obtained from optical and 2MASS NIR data 
respectively. The error bars are derived assuming that the number of stars in each annulus 
follows Poisson statistics. 

The radius of the cluster {rd) is defined as the point where the cluster stellar density 
merges with the field stellar density. The contribution of the field stars has been estimated 
by counting stars in the annular region 16' - 21' from the cluster center. The horizontal 
dashed line in Fig. H] shows the field star density. Optical RDP indicates a dip in stellar 
counts at r ~ 6'. 5, whereas this dip is absent in NIR RDP. The NIR RDP indicates an 
enhancement in stellar density at r ~ 7'. 5. To find out the cause of dip in the optical RDP 
and enhancement in stellar density in the NIR RDP, in Fig. HI we also show RDPs for the 
eastern and western regions of Stock 8. The dip in the optical and enhancement in NIR RDP 
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is visible only towards the western region of the cluster. A comparison of the RDPs with 
spatial distribution oi E{B — V), shown in Fig. [71 indicates that the dip in the optical RDP 
at r ~ 6'. 5 may be an artifact of relatively higher extinction around r ~ 6'. 5 — 7'. The NIR 
observations can probe deeper in the embedded region, hence enhancement in stellar density 
in NIR RDP around r ~ 7'. 5 may be due to embedded population. Since star formation 
activity in the Nebulous Stream, as discussed in Sec. 11, is found to be independent of the 
star formation activity in the cluster region, we used radius of the cluster as ~ 6' (as 
obtained from optical data) to study the properties of the cluster. 

To parametrize the RDP, we fitted the observed RDP with the empirical model of King 
(1962) which is given by. 



where Tc is the core radius at which surface density p(r) becomes half of the central density, 
/o . The best fit to the observed RDPs, obtained by a minimization technique, is shown 
in Fig. m The core radii thus estimated from optical and NIR RDPs are V A ± 0'.2 and 
1'.6 ± O'.l respectively. 

5 REDDENING IN THE CLUSTER 

The extinction in a star cluster arises due to two distinct sources: (1) the general ISM in 
the foreground of the cluster [E{B — V)mm], and (2) the localized cloud associated with the 
cluster [E{B -V) = E{B - V)^ - E{B - V)min], where E{B - V)^ is the reddening of the 
star embedded in the parent cloud. The former component is characterized by the ratio of 
the total-to-selective extinction Ry [= Ay / E{B - V)]= 3.1 (Wegner 1993; He et al. 1995; 
Winkler 1997), whereas for the intracluster regions of young clusters embedded in dust and 
gas cloud the value of Ry may vary significantly (Chini & Wargau 1990; Tapia et al. 1991; 
Hillenbrand et al. 1993; Pandey et al. 2000). 

The ratio of total-to-selective extinction R is an important quantity, that must be known, 
to get the dereddened magnitudes of the stars. Consequently, the value of R affects the 
distance determination significantly. Several studies have already pointed out anomalous 
reddening law with high R values in the vicinity of star forming regions (see. e.g. Neckel & 
Chini 1981, Chini & Krugel 1983, Chini & Wargau 1990, Hillenbrand et al. 1993, Pandey et 
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Figure 4. Stellar density (filled squares with error bars) as a function of radius from the adopted cluster center for the optical 
and 2MASS data. The dashed line represents the mean density level of the field stars and the continuous curve shows the 
least square fit of the King (1962) profile to the observed data points from Kiso Schmidt observations. Open circles connected 
by dotted line and open triangles connected by dot-dashed line represent RDPs for the east and west region of the cluster 
respectively. The error bars represent ± \/~N errors. Box represents the radial range considered as a field region. 



al. 2000, Samal et al. 2007). Since the cluster is embedded in the parent molecular cloud, it 
will be useful to understand the reddening law in and around the Stock 8 region. 

To study the nature of the extinction law in the cluster region, we used two - colour 
diagrams (TCDs) as described by Pandey et al. (2003). The TCDs of the form of (V — A) 
vs. {B — V), where A is one of the broad-band filters {R, I, J, H, K, L), provide an effective 
method for separating the influence of the normal extinction produced by the diffuse inter- 
stellar medium from that of the abnormal extinction arising within regions having a peculiar 
distribution of dust sizes (cf. Chini & Wargau 1990, Pandey et al. 2000). 

The TCDs for the region r < 12' are shown in Fig. [51 The ^j-^^y-j values in the cluster 
region are estimated using the procedure as described in Pandey et al. (2003). The slopes 
of the distributions rriciuster (cf. Pandey et al. 2003) are found to be 1.13 ± 0.01, 1.95 ± 
0.03, 2.41 ± 0.03, 2.56 ± 0.03 for {V - I), {V - J), {V - H), {V - K) vs. {B - V) diagrams 
respectively. Identical values for the slopes are found for the region r < 6'. The ratio of 
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Figure 5. Two - colour diagrams for stars within the radius r 12' from the center of Stock 8. The continuous line represents 
the least square fit to the data. 

total-to-selective extinction in the cluster region, Rdusteri is also derived using the procedure 
given by Pandey et al. (2003). The ratios (A ^ A/) yield Rduster = 3.1 ± 0.1 which 

indicates a normal reddening law in the cluster region. 

In the absence of spectroscopic observations, the interstellar extinction E{B — V) toward 
the cluster region can be estimated using the {U — B)/{B — V) colour-colour (CC) diagram. 
The CC diagram of the cluster region (r ^ 6') is presented in Fig. [6l Since the cluster 
is very young, a variable reddening is expected within the cluster region. The figure also 
indicates a large amount of contamination due to field stars. The probable foreground stars 
follow the intrinsic zero-age-main-sequence (ZAMS) by Girardi et al. (2002) reddened by 
E{B — V)=Q.2Q mag along a normal reddening vector (i.e., E{U — B)/E{B — V) = 0.72). In 
Fig. [6l the continuous line represents the theoretical ZAMS which is shifted hj E{B — V) = 
0.40 and 0.60 mag respectively, along the normal reddening vector to match the observations 
of probable cluster members. Fig. [6] yields a variable reddening with E{B — V)min = 0.40 
mag to E{B — V)max = 0.60 mag in the cluster region. A careful inspection of CC diagram 
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indicates the presence of further reddened population. The theoretical ZAMS, shown by 
dot-dashed line, is further shifted to match the reddened sequence. This population may 
belong to blue plume (BP) of Norma-Cygnus arm (cf. Carraro et al. 2005, Pandey et al. 
2006). The E{B — V) value for the background population comes out to be ~0.70 mag, 
which is comparable to the E{B — V) value of BP population around / ~ 170° (cf. Pandey 
et al. 2006). 

Reddening of individual stars having spectral types earlier than AO have also been com- 
puted by means of the reddening free index Q (Johnson & Morgan 1953). Assuming a normal 
reddening law we can construct a reddening-free parameter index Q = {U — B) — 0.72 x 
E[B — V). For the MS stars, the intrinsic {B — V)q colour and colour-excess can be ob- 
tained from the relation (B — V)q = 0.332 x Q (Johnson 1966, Hillenbrand et al. 1993) and 
E{B — V) = {B — V) — {B — V)o, respectively. The distribution of reddening as a function of 
radial distance from the cluster center is shown in Fig. [TJ which indicates a slight enhance- 
ment of E{B — V^) at r ~ 7'. In Fig. [7] we also show variation of reddening towards eastern 
and western directions of the cluster. From the figure it can be noted that the reddening 
towards eastern region of the cluster is rather uniform whereas the reddening towards the 
western region of the cluster shows relatively higher value (although statistically weak) at 
r ~ 7'. 



6 OPTICAL COLOUR MAGNITUDE DIAGRAMS 

The optical colour- magnitude diagrams (CMDs) were used to derive the cluster fundamental 
parameters such as age, distance etc. The V/ {B — V) and V/ {V — I) CMDs of stars within 
radial distance r ^ 6' of Stock 8 and V/ {V — I) CMD for the nearby field region are shown 
in Fig. [HI The contamination due to background field population is apparent in the CMDs 
of the cluster region. The CMDs of the cluster region also show a significant number of stars 
towards the right of the ZAMS. Stars brighter than V ~ 11.5 mag observed in B-band using 
Kiso Schmidt telescope are saturated even in frames with short exposure. 

To study the LF/MF, it is necessary to remove field star contamination from the sample 
of stars in the cluster region. Membership determination is also crucial for assessing the 
presence of PMS stars because both PMS and dwarf foreground stars occupy similar positions 
above the ZAMS in the CMDs. In the absence of proper motion study, we used statistical 
criterion to estimate the number of probable member stars in the cluster region. To remove 
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Figure 6. {U — B)/{B — V) colour-colour diagram for the stars within the region r ^ 6'. The dotted curve represents the 
ZAMS by Girardi et al. (2002) shifted along the reddening slope of 0.72 (shown as a dashed line) for E{B — V)=0.20 mag 
(probable foreground stars). The continuous curves are the ZAMS shifted by E{B — V) = 0.40 and 0.60 mag respectively to 
match the observations of probable cluster members. The dot-dashed curve represents the ZAMS reddened by E{B — V)=0.70 
mag to match the probable background population (see text). 

contamination of field stars from the MS and PMS sample, we statistically subtracted the 
contribution of field stars from the observed CMD of the cluster region using the following 
procedure. For any star in the V/ {V — I) CMD of the field region, the nearest star in the 
cluster's V/{V-I) CMD within V ±0.125 and (F-/)± 0.065 of the field star was removed. 
While removing stars from the cluster CMD, necessary corrections for incompleteness of the 
data sample were taken into account. The statistically cleaned V/{V — I) CMD of the cluster 
region is shown in Fig. [Hli which clearly shows a sequence towards red side of the MS. The 
contamination due to field stars at V > 20 mag and (V^ — /) ~ 2.1 mag can still be seen in 

Fig. m. 

In Fig. El using E{B- V)m.in = 0.40 mag and relations Ay = 3.1 x E{B -V); E{V - 1) = 
1.25xE{B — V), we visually fitted theoretical isochrone of log age = 6.3 (2 Myr) and Z = 0.02 
by Girardi et al. (2002) to the blue envelope of the observed MS. We estimated a distance 
modulus of (m — M)y = 12.8 ± 0.15 mag which corresponds to a distance of 2.05 ± 0.10 
kpc. Fig. [9] shows dereddened Vq/{B — V)q CMD of the stars lying within r ^ 6'. The stars 
having spectral type earlier than AO were dereddened individually using the Q method as 
discussed in Section |5l Because bright stars in our observations got saturated even in short 
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Figure 7. Distribution of average reddening (filled squares) as a function of radial distance from the cluster center. Error bars 
represent standard errors. Open circles connected by dotted line and open triangles connected by dot-dashed line represent 
distribution of average reddening in the east and west region of the cluster respectively. 



exposure frames, we obtained their B-band data from from Mayer & Macak (1971). The 
isochrone for 2 Myr by Girardi et ah (2002) are also plotted in the figure. This figure yields 
an average post-main-sequence age of the massive stars of the cluster as ^ 2 Myr. 

In Fig. [TU] we have plotted V/ {V — I) CMD for the young stellar objects (YSOs) i.e. Ha 
emission and NIR excess sources (cf. Section 17.11) lying within regions r ^ 6' and 6' < r 
^ 12'. The Ha sources are represented by open and filled triangles and the NIR excess 
sources by open and filled circles. The PMS isochrones by Siess et al. (2000) for 0.5, 5, 10 
Myr (dashed lines) and isochrone for 2 Myr by Girardi et al. 2002 (continuous line) are also 
shown. Fig. [To] reveals that majority of the YSOs (i.e.. Ha and NIR excess stars) have ages 
^ 5 Myr. The age spread may indicate a non-coeval star formation in the cluster. Fig. 
also indicates a tendency that stars lying at r > 6' are relatively younger than those in the 
cluster region. 

In Fig. [11] we have shown statistically cleaned CMD along with PMS isochrones for ages 
0.5 Myr and 10 Myr by Siess et al. (2000). The 2 Myr isochrone by Girardi et al. (2002) and 
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Figure 8. (a) and (b): V/{B - V) and V/{V - I) CMDs for stars within the region r ^ 6' of the cluster Stock 8, (c): V/{V - I) 
CMD for stars in the nearby field region having same area as in Figs.[8ji andfSjj, (d): statistically cleaned V/{V — I) CMD for 
the cluster region. The continuous lino is the isochrone of 2 Myr from Girardi et al. (2002) corrected for the cluster distance 
and reddening. The stars in V/{B — V) CMD shown by open circles are taken from Mayer & Macak (1971). 

evolutionary tracks by Siess et al. (2000) for different masses are also shown. A comparison 
of Figs. [10] and [11] suggests that the stars lying on the right side of the MS in Fig. [11] may 
be probable PMS stars. The stars lying between the 0.5 - 10 Myr isochrones are considered 
as PMS stars and this set of data is used for calculation of the IMF (see Sec. 8). 

7 NEAR INFRARED COLOUR-COLOUR AND COLOUR-MAGNITUDE 
DIAGRAMS 

The NIR data are very useful to study the nature of young stellar population within the 
star forming regions. From 2MASS point source catalogue, we obtained JHKg data (with 
photometric errors ^0.1 mag in all the three bands) for 445 and 898 NIR sources lying within 
r ^ 6' and 6' < r ^ 12' respectively. In the following Section, we discuss NIR colour-colour 
diagram and CMDs. 
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Figure 9. Vo/{B — V)q CMD for stars lying within the region r ^ 6' of the cluster Stock 8. The data for stars shown by 
open circles are taken from Mayer & Macak (1971). The isochrone (continuous curve) for 2 Myr by Girardi et al. (2002) is also 
shown. The average photometric errors in magnitude and colour are shown in the lower left of the figure. 

7.1 Colour-Colour Diagram 

The (J — H)/{H — K) colour-colour diagrams for the cluster region r ^ 6', outside the 
cluster region 6' < r ^ 12', and a nearby field region are shown in Fig. [121 The thin and 
thick solid lines are the locations of unreddened main-sequence and giant stars (Bessel & 
Brett 1988), respectively. The dotted line represents the locus of T Tauri stars (Meyer et al. 
1997). The parallel dashed lines are the reddening vectors for the early MS and giant type 
stars (drawn from the base and tip of the two branches). The crosses on the dashed lines are 
separated by an Ay value of 5 mag. The extinction ratios, Aj/Ay = 0.265, Ah/ Ay = 0.155 
and Ax/Ay = 0.090, are adopted from Cohen et al. (1981). All the 2MASS magnitudes and 
colours have been converted into the CIT system. The curves are also in the CIT system. 

Presently young stellar objects are classified as an evolutionary sequence spanning a 
few million years as: Class 0/ Class 1 - the youngest embedded protostars surrounded by 
infalling envelopes and growing accretion disk; Class 11 - PMS stars with less active accretion 
disks and Class III - PMS stars with no disks or optically thin remnant dust (Adams et al. 
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Figure 10. The V/{V — I) CMD for Ha emission (open and filled triangles) and NIR excess stars (open and filled circles), 
lying within the regions r 6' and 6' < r ^ 12' respectively, fsochrone for 2 Myr by Girardi et al. (2002) (continuous curve) 
and PMS isochrones for 0.5, 5, 10 Myr by Siess et al. (2000) (dashed curves) are also shown. All the isochrones are corrected 
for the distance of 2.05 kpc and reddening E(B — V) = 0.40 mag. Average photometric errors in {V — I) colour as a function 
of magnitudes are shown in the left side of the figure. 



1987). Following Ojha et al. (2004a), we classified sources according to their locations in 
(J — H)/{H — K) colour-colour diagrams. The 'F' sources are those located between the 
reddening vectors projected from the intrinsic colour of main-sequence and giant. These are 
considered to be field stars (main-sequence, giants) or Class III / Class II sources with small 
NIR excesses. The 'T' sources are located redward of 'F' but blueward of the reddening line 
projected from the red end of the T Tauri locus of Meyer et al. (1997). These sources are 
considered to be mostly classical T Tauri stars (Class II objects) with large NIR excesses. 
There may be an overlap in NIR colours of Herbig Ae/Be stars and T Tauri stars in the 'T' 
region (Hillenbrand et al. 1992). The 'P' sources are those located in the region redward of 
region 'T' and are most likely Class I objects (protostellar objects). 

In Fig. [121 ^ significant number of stars show NIR excess indicating that the Stock 8 
region is populated by YSOs. A comparison of the colour-colour diagrams of the cluster 
region with the field region (Fig. [T2b ) indicates that stars in the cluster region having 
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Figure 11. Statistically cleaned V/ (V — I) CMD for stars lying within the region r 6' of Stock 8 and between PMS isochrones 
of age 0.5 Myr and 10 Myr. The isochrone for 2 Myr age by Girardi et al. (2002) and PMS isochrones of 0.5, 10 Myr along with 
evolutionary tracks of different mass stars by Siess et al. (2000) are also shown. All the isochrones are corrected for the cluster 
distance and reddening. The corresponding values of masses in solar mass are given at the right side of each track. Points shown 
by small dots are considered as non-members. Average photometric errors in {V — I) colour as a function of magnitudes are 
shown in the left side of the figure. 



J — H > 0.6 and lying towards the right side of the reddening vector at the boundary of 
'F' and 'T' regions can be safely considered as NIR excess stars. This criterion yields 19 
and 17 NIR excess sources in the Stock 8 (r ^ 6') and outside the cluster (6' < r ^ 12') 
regions, respectively, which are shown as open and filled circles in Fig. [121 Distribution of 
Ha emission stars is also shown as open (r ^ 6') and filled (6' < r ^ 12') triangles. Fig. [12] 
indicates that the YSOs located outside the cluster region have relatively higher extinction 
and NIR excess as compared to those located in the inner region. The distribution of NIR 
excess, A(i/ — K), defined as horizontal displacement from the reddening vector at the 
boundary of 'F' and 'T' regions (cf. Fig. [T2|) . for the two regions is shown in Fig. [T3] which 
manifests that YSOs located in the region 6' < r ^ 12' have relatively higher NIR excess. 
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Figure 12. {J — H)/{H ~ K) colour-colour diagrams of sources detected in the JHKs bands with photometric errors less than 
0.1 mag in (a) Stock 8 region (r 6') (b) outside the Stock 8 region 6' < r ^ 12' (c) the nearby reference field. The locus for 
dwarfs (thin solid curve) and giants (thick solid curve) are from Bessell & Brett (1988). The dotted line represents the locus 
of T Tauri stars (Meyer et al. 1997). Dashed straight lines represent the reddening vectors (Cohen et al. 1981). The crosses on 
the dashed lines are separated by Ay = 5 mag. The open and filled triangles are Ha emission sources while open and filled 
circles are NIR excess sources for the inner and outer regions, respectively. The average photometric errors arc shown in the 
lower left of the panel. 



7.2 Colour-Magnitude Diagrams 

The K/ {H — K) and J / (J — H) CMDs for the NIR excess and Ha sources detected within 
r ^ 12' region are shown in Figs. [T3] and [151 The isochrone for age 2 Myr by Girardi et al. 
(2002) and PMS isochrones for ages 0.5, 10 Myr by Siess et al. (2000) respectively, have been 
plotted assuming a distance of 2.05 kpc and an extinction of E{B — V)min = 0.40 mag as 
obtained from the optical data. Fig. [T3] also indicates that sources lying outside the cluster 
region i.e. 6' < r ^ 12' have relatively higher NIR excess. Therefore, from Figs. [TUl [T21 [13] 
and [131 we infer that sources lying in the region 6' < r ^ 12' may possibly be younger than 
those located in the cluster region (r ^ 6'). 

The mass of the probable YSO candidates can be estimated by comparing their locations 
on the CMD with the evolutionary models of PMS stars. To estimate the stellar masses, the 
J luminosity is recommended rather than that of H or as the J band is less affected by 
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Figure 13. Distribution of NIR excess in two regions. Continuous and dashed line histograms show the distribution in r 6' 
and 6' < r ^ 12' regions respectively. 



the emission from circumstellar material (Bertout et al. 1988). In Fig. [151 the continuous 
obhque reddening hues denote the positions of PMS stars of 0.5 Myr age having masses 0.1, 
1.0 and 3.0 Mq. The YSOs, in general, have masses in the range 0.1 to 3.0 Mq. 



8 INITIAL MASS FUNCTION 

Young clusters are important tools to study the IMF since their MF can be considered as 
the IMF as they are too young to loose significant number of members either by dynamical 
or stellar evolution. To study the IMF of Stock 8 we used the data within r ^ 6'. 

The mass function (MF) is often expressed by the power law, A^(logm) oc mF and the 
slope of the MF is given as: 

r = (ilogA'"(logm)/(ilogm 

where A^(logm) is the number of stars per unit logarithmic mass interval. For the mass 
range 0.4 < M/Mq ^ 10, the classical value derived by Salpeter (1955) for the slope of IMF 
is F = -1.35 . 
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Figure 14. K/{H — K) CMD for stars within the region r ^ 12'. The symbols are same as in Fig. 1101 Isochrone for 2 Myr 
by Girardi ct al. (2002) (continuous curve) and PMS isochrones from 0.5 and 10 Myr by Siess et al. (2000) (dashed curves), 
corrected for cluster distance and reddening are also shown. 



With the help of statistically cleaned CMD shown in Fig. [TT|, we can derive the MF 
using theoretical evolutionary models. Since the post-main-sequence age of the cluster is 
^ 2 Myr, stars having ^ 15 mag have been considered to be on the main sequence. 
For the MS stars LF was converted to MF using the theoretical model by Girardi et al. 
(2002) (cf. Pandey et al. 2001, 2005). The MF for PMS stars was obtained by counting the 
number of stars in various mass bins (shown as evolutionary tracks in Fig. [TT]) . Necessary 
corrections for incompleteness of the data sample were taken into account for each magnitude 
bin to calculate the MF. The MF of the cluster is plotted in Fig. (TBI The slope of the mass 
function F, in the mass range 1.0 ^ M/Mq < 13.4, can be represented by a power law 
having a slope of F = —1.38 ± 0.12, which agrees well with Salpeter value (-1.35). In the 
mass range 0.3 ^ M/Mq < 1.0, the mass function also follows a power law but with a 
shallower slope F = —0.58 ± 0.23 indicating a break in the slope of the MF at ~ IMq. 
The break in the power law has already been reported in the case of a few young clusters 
e.g. Trapezium and IC 348 (Muench et al. 2002, 2003). They showed a flattening in slope 
of the IMF at ~ 0.6Mq. Recently, in the case of young clusters NGC 1893 (Sharma et al. 
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Figure 15. J/{J — H) CMD for stars within the region r 5C 12' . The symbols are same as in Fig. [To] The isochrone of 2 Myr 
(continuous curve) by Girardi et al. (2002) and PMS isochrones of age 0.5 and 10 Myr (dashed curves) by Siess et al. (2000), 
corrected for cluster distance and reddening arc also shown. The continuous oblique reddening lines denote the positions of 
PMS stars of 0.1, 1.0 and 3.0 A/q for 0.5 Myr. 

2007) and Be 59 (Pandey et al. 2007), a break in the power law was reported at ~ 2Mq 
and ~ 2.5Mq, respectively. However, F = —0.58 ± 0.23 obtained for Stock 8 below IMq is 
shallower than that obtained for NGC 1893 (—0.88 ±0.09) and is steeper than that obtained 
for Be 59 (approximately flat). In the case of NGC 2362 the MF flattens at ~ 3 - IMq 
(Damiani et al. 2006) which is at a much higher mass than observed in some young clusters. 
Prinsinzano et al. (2005) found that the MF of young open cluster NGC 6530, in the mass 
range 0.4 ^ M/Mq < 4.0, can be represented by a power law having slope (—1.22 ± 0.17). 
However the MF does not flatten below ~ O.6M0 as in the case of some young clusters e.g. 
Trapezium, Taurus and IC 348, instead it decreases for masses lower than ~ O.4M0. Here 
it is worthwhile to mention that exclusion of the lowest mass bin in Fig. [16] will yield a 
shallower slope F = —0.3 ± 0.3, although with large error, for masses below IM0, which is 
comparable with the slope obtained by Muench et al. (2002, 2003) for Trapezium and IC 
348 clusters. However a deeper photometry is required to study the nature of mass function 
of Stock 8 below IM^. 
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Figure 16. A plot of the mass function for the cluster region r ^ 6'. The log (j> represents A^/dlog m. The error bars represent 
errors. The continuous lines show least-squares fit to the mass ranges described in the text. The values of the slopes 
obtained are also mentioned in the figure. 

The mass segregation, in the sense that massive stars being more centrally concentrated 
than the lower mass stars, has been reported in several LMC and Milky Way star clusters 
(see Chen et al. 2007; de Grijs et al. 2002 a, b, c; Fischer et al. 1998; Pandey et al. 1992, 
2001, 2005). In the case of intermediate or old clusters, the mass segregation is mainly due 
to the effect of dynamical evolution. But in the case of young clusters the mass segregation 
may be the imprint of star formation process. To characterize the degree of mass segregation 
in Stock 8, we plotted cumulative distribution of stars as a function of distance from the 
cluster center in two mass groups as shown in Fig. [T71 The left panel shows the distribution 
for mass groups 4 ^ M/Mq < 13.4 (high mass group) and 0.3 ^ M/Mq < 4 (low mass 
group), whereas the right panel shows the distribution for 1.0 ^ M/Mq < 13.4 (high mass 
group) and 0.3 ^ M/Mq < 1.0 (low mass group). The Kolmogorov-Smirnov test indicates 
that the left panel shows mass segregation at a confidence level of ~ 75%, whereas the right 
panel does not show any mass segregation. This indicates that mass segregation, although 
weak, is effective only towards the higher mass end. Since variation of reddening within the 
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Figure 17. Cumulative distributions of stars of different mass intervals as a function of radial distance. 



cluster is not significant (cf. Fig. [7]), and incompleteness of the data does not depend on the 
radial distance (cf. Table [3]), we believe that the mass segregation may be the imprint of 
star formation process. 

The i^-band luminosity function (KLF) is a powerful tool to investigate the IMF of 
young embedded star clusters. During the last decade several studies have been carried out 
with the aim of determining the KLF of young open clusters (e.g. Lada & Lada 2003, Ojha et 
al. 2004b, Sanchawala et al. 2007). We have used 2MASS Kg-band data to study the KLF in 
Stock 8 region. The completeness of the 2MASS data is obtained by using ADDSTAR routine 
as discussed in Sec. 2.1. The CFs obtained for two sub-regions are given in Table [51 The 
KLF in Stock 8 region is studied by using the Besangon Galactic model of stellar population 
synthesis (Robin et al. 2003) and stars from a nearby reference field to take into account 
foreground/background field star contamination (cf. Sharma et al. 2007). The reference field, 
also taken from 2MASS, is located at a radial distance of ~ 30' away from the cluster center. 
Star counts are predicted using the Besangon model towards the direction of the control field. 
An advantage of using this model is that we can simulate foreground {d < 2.05 kpc) and the 
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Figure 18. KLF derived after subtracting the model corrected field star contamination (see the text). The linear fit is 
represented by the continuous line. 



background {d > 2.05 kpc) field star populations separately. The foreground population was 
simulated using the model with Ay = 1.24 mag {E{B — V^) = 0.40 mag (cf. Section [5]) and 
d < 2.05 kpc. The background population {d > 2.05 kpc) was simulated with Ay = 1.86 
mag. Thus we determined the fraction of the contaminating stars (foreground + background) 
over the total model counts. This fraction was used to scale the nearby reference field and 
subsequently the star counts of the modified control field were subtracted from the KLF of 
the cluster to obtain the final corrected KLF. The scale factor in different magnitude bins 
was found to be in the range of 0.8 - 1.0 with a mean value of ~ 0.9. 
The KLF is expressed by following power-law: 



dN(K) 



oc 10 



aK 



dK 

where is the number of stars per 0.5 magnitude bin and a is the slope of the power 
law. The KLF for the cluster region r ^ 6' shown in Fig. [181 yields a slope of 0.31 ± 0.02, 
which is smaller than the average value of slopes {a ~ 0.4) for young clusters (Lada et al. 
1991; Lada & Lada 1995; Lada & Lada 2003). Smaller values of KLF slope (~ 0.3 — 0.2) have 
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Table 5. Completeness Factor (CF) of 2MASS data in the cluster and field regions. 



K range 


Stocks 




Field region 


(mag) 


r ^ 3' 


3' < r ^ 6' 




10.0 - 10.5 


1.00 


1.00 


1.00 


10.5 - 11.0 


1.00 


1.00 


1.00 


11.0 - 11.5 


1.00 


1.00 


0.98 


11.5 - 12.0 


1.00 


1.00 


0.97 


12.0 - 12.5 


0.95 


0.96 


0.91 


12.5 - 13.0 


0.99 


1.00 


0.96 


13.0 - 13.5 


0.94 


0.94 


0.97 


13.5 - 14.0 


0.96 


0.96 


0.97 



been reported for various young embedded clusters (Megeath et al. 1996, Chen et al. 1997, 
Brandl et al. 1999, Ojha et al. 2004b, Leistra et al. 2005, Sanchawala et al. 2007, Pandey et 
al. 2007). 

The KLF of Stock 8 is worth comparing with those of NGC 1893 and Be 59 since all the 
KLFs are obtained using a similar technique (NGC 1893, Sharma et al. 2007; Be59, Pandey 
et al. 2007). The slope {a = 0.31 ± 0.02) obtained for Stock 8 is in agreement with those 
obtained for NGC 1893 (a = 0.34 ± 0.07) and Be 59 (a = 0.27 ± 0.02). 

9 EMBEDDED CLUSTER IN THE NEBULOUS STREAM 

A Nebulous Stream, as identified in Fig. [H is clearly seen towards the eastern region of the 
cluster. Stock 8. An embedded cluster, at a radial distance of ~ 13' from the center of Stock 
8, can be noticed in the Nebulous Stream at aaooo = 05^^29^00^; ^aooo = +34°23'12". Its 
2MASS Kg band image together with the distribution of NIR excess stars and Ha emitters 
is shown in Fig. [191 Borissova et al. (2003) detected this cluster in their search for 2MASS 
PSC and designated it as CC 14. Recently, Ivanov et al. (2005) have studied CC 14 using 
deep NIR observations. They estimated a distance modulus (m — M)o = 16.6 ± 0.8 mag, 
corresponding to a distance of ~ 21 kpc which indicates that CC 14 is not associated with 
the IC 417 region. However, the distribution of NIR-excess and Ha emission stars shown in 
Fig. [19] seems to indicate a physical relationship between CC 14 and the H II region IC 417, 
which is discussed later. In Fig. [20] we show the RDP for CC 14 obtained from the 2MASS 
data around the position mentioned above, which indeed manifests a clustering with a radius 
of ~ 1'.5. 

In Fig. [2Tk we show (J — H) / [H — K) colour-colour diagram for stars in CC 14 cluster as 
well as for the YSOs distributed along the Nebulous Stream. A comparison of Fig. [2Tk with 
Fig. [12] indicates that the stars having (J — H) > 0.6 can be considered as probable YSOs 



30 Jessy Jose et al. 




Figure 19. Spatial distributions of Ha emitters (open circles) and IR-excess sources (asterisk symbols) overlaid on DSS-2 R 
band image. The thick contours represent ^^CO emission from Leisawitz et al. (1989). NVSS (1.4 GHz) radio contours (white) 
and MSX A-band intensity contours (thin black) have also been shown. The MSX A-band contours are 1.5, 2, 3, 4, 5, 10, 20, 
40, 60, 80 % of the peak value 9.666 X 10"'' Wm.-'^Sr-^ and the NVSS radio contours are 3, 5, 10, 20, 40, 60, 80, 90 % of the 
peak value 29.5 mjy/beam. The circle having radius, r = 6' represents optical extent of the cluster (see Fig.|4j. The inset box 
shows the enlarged 2MASS A's-band image of the embedded cluster CC14. The abscissa and the ordinates are for the J2000 
epoch. 



associated with CC 14 and the Nebulous Stream. Assuming that the YSOs have average 
mass of ~ (i.e. My ~ 4.5 at 1 Myr), mean V ~ 18.5 mag and average Ay ~ 2.5 mag 
(cf. Fig. [2T!) . we estimate a true distance modulus of ~ 11.5, which clearly indicates that 
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Figure 20. Radial density profile for the embedded cluster CC14 obtained from the 2MASS data. The continuous curve is the 
least square fit of the King (1962) profile to the observed data points and the error bars represent it \/N errors. The dashed 
line represents the mean density level of the field stars. 



CC 14 and YSOs associated with the Nebulous Stream can not be located at the distance 
(~ 21 kpc) suggested by Ivanov et al. (2005). 

Of the 23 stars located within CC 14 region, the location of 10 stars (open circles) on 
the (J — H) / [H — K) colour-colour diagram indicates that they may be probable field stars. 
The RDP (Fig. [20i) for CC 14 also indicates the presence of ~ 14 field stars in the region. 
Probable members of CC 14 are shown by filled circles in Fig. I^Tk . Figures [?Tb and [?Ib 
show V/{y — I) and J/ {J — H) CMDs for the CC 14 region where open and filled circles 
represent field and probable cluster populations respectively. Comparison of V/{V — I) CMD 
(Fig. Elb) with the field region V/{V - I) CMD (Fig. [Hb) clearly indicates that the field 
population (open circles) in the embedded cluster region belongs to the same population 
as that of the nearby field region. The probable cluster members (filled circles) follow ~ 1 
Myr isochrone. In Figs. [2Tk-c we have also plotted the NIR-excess and Ha emission stars 
lying along the Nebulous Stream using crosses and star symbols, respectively. The V/{V — I) 
CMD indicates that the stars lying on the Nebulous Stream also have ages of ~ 1 Myr. 
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Figure 21. Colour - colour diagram and CMDs for the stars lying within 1'.5 radius of the embedded cluster CC14 as well 
as for YSOs lying along the Nebulous Stream, (a): (J — H)/{H — K) colour - colour diagram with probable cluster members 
(J — H > 0.6) as filled circles and field stars as open circles. Crosses and star symbols are the NIR excess and Ha emission 
sources respectively, lying along the Nebulous Stream. (6) and (c) are V/{V — I) and J/{J — H) CMDs respectively. The symbols 
are same as in Fig. 121b .. In Figs. 121b and I21t the isochrone of 2 Myr (continuous curve) by Girardi et al. (2002) and PMS 
isochrones of age 1 and 10 Myr (dashed curves) by Siess et al. (2000), corrected for cluster distance and reddening, are also 
shown. 

10 DISTRIBUTION OF GAS AND DUST AROUND STOCK 8 

No detailed observations have been made so far, to probe the molecular cloud associated 
with Stock 8. Only available is the ^^CO contour map of IC 417 (Sh2-234) region traced 
from figure 27c of Leisawitz et al. (1989), which is shown in Fig. [191 Although the resolution 
is low, a moderate-sized molecular cloud of ~ 3.4 x IO^Mq is found adjacent to the east of 
Stock 8. 

In Fig. [191 NVSS (1.4 GHz) radio contours (white) and MSX A-band MIR contours 
(black) are overlaid on DSS-2 R band image. The radio continuum emission peaks around 
«2ooo ~ 05^*28'" 18" and ^2000 ~ +34°25'30" and shows a sharp cutoff at aaooo ~ 05''28'"30" 
indicating an interface between the ionized gas and molecular cloud. The integrated flux 
density of the radio continuum above 3a level is estimated to be 2.474 Jy. Assuming a 
spherical symmetry for the ionized region and neglecting absorption of ultraviolet radiation 
by dust inside the H II region, the above flux density together with assumed distance, allow 
us to estimate the number of Lyman continuum photons {Niyc) emitted per second, and 
hence the spectral type of the exciting star. Using the relation given by Martin-Hernandez et 
al. (2003) and assuming an electron temperature of 8000 K, we estimated logN^yc = 47.67, 
which corresponds to a MS spectral type of ~ BO (Panagia 1973; Schaerer & de Koter 1997). 
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Figure 22. Spatial distributions of OB stars, Ha emitters (open circles) and IR-excess sources (asterisks) overlaid on DSS-2 
R band image of Stock 8 region. The inset box shows an enlarged view of 6' central region of Stock 8. The thick contours 
represent ^-^CO emission from Leisawitz et al. (1989). MSX A-band intensity contours (thin lines) have also been shown. The 
contour levels arc the same as in Fig. 1191 The star numbers from Table 3 and three sub regions (I - III, see the text) are also 
shown. 

This seems to be consistent with the BO/09.5 star (cf. Table 3, star no. 1), at the center 
of Stock 8, as the ionizing source of the region. Major uncertainty in the above calculation 
may arise due to neglect of the absorption of Lyman continuum photons by the dust. We 
estimated the effect of dust using the following two methods; 1) the dust continuum density 
was estimated by averaging FIR optical depth rioo (shown in Fig. [25h) over the radio emission 
region, and 2) attributing the spread in E{B — V) in the cluster region (cf. Sec. 5) to the 
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Table 6. List of OB stars in IC 417 region 



Star 


0^(2000) 


'^(2000) 


Spectral 


E{B — V) 


References 




no. 


(h:m:s) 


{d:m:s) 


type 


mag 








1 


05:28:07 


+34:25:27 


09.5 V 




GeorgcL 


in, et al.. 


1973 








BO IV 


0.48 


Mayer ^ 


z Macak, 


1971 


2 


05:27:35 


+34:27:00 


08 V 


1.25 


Mayer & 


z Macak, 


1971 


3 


05:27:29 


+34:25:03 


09 IV-V 


0.98 


Mayer & 


z Macak, 


1971 


4 


05:27:23 


+34:23:41 


BO V 


0.55 


Mayer ^ 


z Macak, 


1971 


5 


05:28:08 


+34:25:14 


Bl V 


0.46 


Mayer & 


z Macak, 


1971 


6 


05:27:06 


+34:19:32 


BO. 5 V 


0.44 


Mayer & 


z Macak, 


1971 


7 


05:27:16 


+34:25:44 


BO. 5 V 


0.45 


Mayer ^ 


z Macak, 


1971 


8 


05:28:08 


+34:23:45 


BO. 5 V 


0.47 


Mayer & 


z Macak, 


1971 


9 


05:28:06 


+34:27:22 


B1.5 V 


0.52 


Mayer ^ 


z Macak, 


1971 


10 


05:27:17 


+34:30:57 


B2 V 


0.54 


Mayer ^ 


z Macak, 


1971 


11 


05:27:43 


+34:31:56 


BO. 5 IV 


0.60 


Savage, 


1985 




12 


05:28:06 


+34:25:00 


Bl V 


0.45 


Mayer ^ 


z Macak, 


1971 


13 


05:28:39 


+34:40:09 


08 V 


0.57 


Savage, 


1985 





local dust distribution. In both the cases, we assumed a uniform dust distribution and the 
dust composition to be an equal mixture of Silicate and Graphite with properties identical 
to those given by Draine & Lee (1984). The calculation yields logA'^^^j^c ~ 48.15 (from method 
1) and 48.32 (from method 2) implying the ZAMS spectral type of the ionizing star to be 
09/08.5 and 08 from methods 1 and 2, respectively. Keeping in mind the uncertainties in 
estimating the spectral type of the star and the assumptions made to calculate the total 
number of Lyman continuum flux, we feel that the above estimates are comparable. 

However star no. 1 is located almost at the western edge of the ionized gas, so roughly 
half of its UV emission may not be used in creating the HII region. The H II region IC 417 
contains several OB stars (cf. Mayer & Macak, 1971), whose distribution is shown in Fig. 
I22l Presumably these stars are the first generation of stars in the IC 417 region. The details 
of the OB stars are given in Table 3. We suspect that the 08 and 09 stars (star Nos. 2 
and 3, respectively), which jointly emit Lyman continuum photons of log Niyc = 48.78, also 
contribute to ionize the IC 417 HII region to a similar extent of that by BO/09.5 star in 
Stock 8, although they are located relatively far away. 

The emission from MSX A-band is completely absent from the inner region (r < 2') 
as well as from the western region of Stock 8. It is stronger and more extended towards 
the eastern side of the cluster and then protrudes further to the east, nicely following the 
Nebulous Stream and showing a slight enhancement around CC 14. This further supports 
the physical connection between the H II region and the Nebulous Stream/CC 14. We show 
the MSX A, C, D and E-band contour maps of the Stock 8 region in Fig. [231 The contours 
of the MSX A and C-bands show an elliptical ridge around ^2000 ~ 05^^28^26*, ^2000 ~ 
+34°25'00". The eastward contours seen in A (8 fim) and C (12 fim) bands lie beyond 
the ionization front and get diluted as we move towards D (15 fim) and E (21 fim) bands. 
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Figure 23. (top left): MSX A-band intensity contours with levels as 80, 60, 40, 20, 10, 5, 4, 3, 2, 1.5 % of the peak value 
9.666 X IQ-^Wm-^Sr-''-. {top right) and (bottom left): C and D-band intensity contours with levels as 90, 80, 60, 40, 20, 15, 
12 % of the peak value 1.218 X IQ-^Wm-^Sr-^ and 7.239 X IQ-'^Wm-'^ Sr'^ respectively, (bottom right): E-band intensity 
contours with levels as 90, 80, 70, 60, 50, 40, 35 % of the peak value 6.298 X 10~^Wm~-^ Sr"^ . The abscissa and ordinates are 
in J2000 epoch. 

whereas the emission towards the west remains present in all the four bands. The emission in 
the MIR bands traces presence of warm small dust grains. The bands A and C also include 
several PAH features as described in Sec. 3. The morphology of the emission shown in Fig. 
123) suggests that the westward emission may be due to warm small dust grains, whereas 
eastward emission in A and C bands may be due to PAHs because PAHs cannot survive in 
ionized regions as the hard radiation field within HII regions destroys them (Cesarsky et al. 
1996). 
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Figure 24. IRAS-HIRES intensity maps of 12 fim {top left), 25 fim [top right), 60 fim [bottom left) and 100 (bottom 
right). The contours are at 3, 4, 5, 7, 10, 20, 40, 60, 80, 90 % of the peak value of 190.5 MJy/Sr at 12 fim; 7, 10, 20, 30, 40, 
50, 60, 70, 80, 90 % of the peak value of 97.7 MJy/Sr at 25 fim; 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, 95 % of the peak values 
of 332.6 MJy/Sr and 515.4 MJy/Sr at 60 fim and 100 fim respectively. The abscissa and ordinates are in J2000 epoch. 

In Fig. [21] we show IRAS-HIRES intensity maps for the cluster region at 12/im, 25/im, 
60/xm and lOOyum. The emission at FIR (60/im and 100/xm) peaks at 02000 ~ 05'^28™24''; 
^2000 ~ +34°24'00", whereas at MIR (25/im) the emission peaks at ^2000 ~ 05'^28'"10''; 
<^2ooo ~ +34°25'30". The peaks of FIR and MIR emissions are ~ 2' away from the peak 
of the radio continuum emission, which indicates an interaction between the gas and dust 
around the region. IRAS-HIRES maps at 60 and lOO/im were also used to generate the 
spatial distribution of dust colour temperature (T(60/100)) and optical depth maps at 100 
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Figure 25. {left): The dust optical depth (rioo at 100 /im ) distribution of the region around the cluster. The contours are at 
30, 40, 50, 60, 70, 80, 90, 95 % of the peak value 6.85 X lO"". {right): Dust colour temperature map of the region around the 
cluster. The contour levels are at 34, 36, 39, 42, 44, 48 K. 

yum (tioo) using the procedure given by Ghosh et al. (1993). An emissivity law of ex cx A"^ was 
assumed to generate the optical depth map. Fig. [25] presents the contour maps of dust colour 
temperature (T(60/100)) and optical depth (rioo)- The dust colour temperature varies from 
34 to 48 K as can be inferred from the map shown in Fig. [251 The relatively high temperature 
around the cluster may probably be due to the radiation from massive star(s) in the region. 
The Tioo map represents low optical depth around the position of Stock 8, whereas its peak 
matches well with the CO peak and the position of CC 14. 

11 NATURE OF THE NEBULOUS STREAM 

What is the nature of the Nebulous Stream? At first impression it appears like a giant 
Herbig-Haro jet and reminds us of HH 399 in Trifid Nebula (Yusef-Zadeh et al. 2005). 
However this does not seem to be the case. Fig. [19] reveals that it is closely associated with 
the distribution of the MSX A-band emission, as well as the distribution of Ha emission 
stars and IR-excess stars. The cluster CC 14 is also embedded within the Nebulous Stream. 
We therefore interpret the Nebulous Stream as a bright rim or an ionization front. 

Fig. [26] shows continuum-subtracted Ha and [S II] images for the region along with 



38 Jessy Jose et al. 




Figure 26. Continuum-subtracted (20 X 10 arcmin^) Ha {upper panel) and [S II] (lower panel) images of the regions around 
Stock 8 overlaid with MSX A- band contours (thick black) and NVSS radio continuum contours (Ught black). The contour 
intensity levels are the same as mentioned in Fig. 1191 The abscissa and the ordinates are for the J2000 epoch. The discontinuity 
in the images at 02000 ~ 05''28'"30'' is an artifact of observations as the region shown in the figure is mosaic of two non- 
overlapping frames. The area between dashed lines represent a 0.5 arc-min wide strip along the Nebulous Stream (see text). 
Spatial variation of [S II] /Ha ratio along the strip is shown in Fig. 1271 



contours for MSX A-band emission and radio continuum superimposed on the [S II] image. 
The gap seen at the middle of Fig. [26] is because the field was observed in two frames without 
an overlap. The Ha emission is more extended compared to the [S II] emission. The [S II] 
emission comes from low-excitation zones and is enhanced close to the ionization fronts (cf. 
Deharveng et al. 2003). Therefore the spatial distribution of the [S II] /Ha ratio should give 
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Figure 27. Spatial variation of [S II] /Ha ratio within a 0'.5 wide strip along the Nebulous Stream from Stock 8 as shown in 
Fig. 1261 The error bars represent the standard deviation. 



information on the excitation conditions. The variation of the [S II] /Ha ratio along a 0'.5 
wide strip (indicated by black dashed lines in Fig. [26l) as a function of radial distance from 
Stock 8 is shown in Fig. [271 To improve the signal-to-noise ratio and to remove faint stars, 
the images were median filtered using a 3 x 3 pixel box. A background was then determined 
for these smoothed images by calculating the statistics inside a 3x3 pixel box at various 
places where the Ha and [S II] brightness is at a minimum. The average background value 
was then subtracted from each pixel. As can be seen in Fig. [27] the ratio of [S II] /Ha is 
higher eastwards of the cluster Stock 8. Note that the values in the eastern most region are 
not reliable because the error bars are very large. Mean value of the [S II] /Ha ratio for the 
central region r < 6' of Stock 8 is ~ 0.18±0.06, which is consistent with the ratio expected 
for an H II region associated with a late 0-BO type star (Reynolds 1988). Whereas in the 
region between 7' and 12 ' the mean value of the ratio is higher, ~ 0.24±0.02. This indicates 
that the Nebulous Stream is a low-excitation object suggestive of its bright-rim/ ionization 
front nature. 

The question then arises; where is(are) the source(s) of its excitation? The 09.5/ BO 
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stars of Stock 8 or 08/ 09 stars further west do not seem to be the possible source because 
the morphology of the radio and MSX A-band emissions indicates an interface at 02000 ~ 
05''28.5™ and suggests that the ionization front has not reached the Nebulous Stream. UV 
radiation from the west of the Stream is most probably blocked by the molecular gas located 
at the interface at ^2000 ~ 05^*28. 5"\ In addition, the Stream runs in the east-west direction, 
which suggests that the UV photons responsible for the Stream should arrive from the north. 
In a search towards the north of the Stream, we find an 08 star (star number 13 at 0:2000 = 
05''28™39^ (52000 = +34°40'09", distance ~ 2.06 kpc: Cruz-Gonzalez et al. 1974; 2.7 kpc: 
Savage, 1985) located at a projected distance of 9 pc from the Nebulous Stream. We suspect 
that this star could be the possible source of UV photons which excite the Stream. 

In order to verify whether this star emits enough UV photons to generate the bright 
rim, we calculated the flux of Lyman continuum photons which are expected to arrive at 
the location of the bright rim. The flux associated with an 08 star is taken from Panagia 
(1973). We assumed in our calculations, that loss of photons due to absorption by interstellar 
matter between the star and bright rim is negligible. Using a projected distance of ~ 9 pc 
between the star and the bright rim, the Lyman photon flux at the bright rim is estimated 
to be ~ 4 X 10^ cm~^ s^^. The Lyman photon flux reaching a radial distance of 6'(~ 3.6 pc) 
from BO star is ~ 2.7 x 10^ cm~^ s~^. The estimated flux of ~ 4 x 10*^ cm~^ in the case 
of 08 star is comparable to the ionizing flux associated with the bright rimmed clouds (see 
Thompson et al. 2004, Morgan et al. 2004) which supports an 08 star as responsible for the 
Nebulous Stream. 

12 STAR FORMATION SCENARIO IN THE VICINITY OF STOCK 8 

The energetic stellar UV radiation and winds from massive stars could disperse nearby clouds 
and consequently terminate further star formation. Alternatively, they can enhance star 
formation activity in the region. Elmegreen & Lada (1977) proposed that after formation of 
massive stars, the expanding ionization fronts play a constructive role in inducing a sequence 
of star-formation activities in the neighbourhood. The distribution of visible young stars 
and embedded YSOs and the morphological details of the environment around the cluster 
containing OB stars can be used to probe the star formation history of the region. 

In addition to the ^^CO contours and the locations of OB stars. Fig. [22] also shows 
spatial distribution of NIR-excess sources and Ho emission stars. In Fig. [28] we plot NIR 
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Figure 28. Variation of NIR excess {A{H — K)) as a function of radial distance from the center of Stock 8 for three sub-regions 
as shown in Fig. 1221 The triangles, circles and asterisks are the sources falling in regions I, II and III respectively. 



excess, A{H — K), for three sub-regions (regions I, II and III) shown in Fig. |22]as a function 
of radial distance from the center of the cluster Stock 8. Fig. [2H] indicates that the stars 
lying in regions I and III show relatively higher NIR excess. This trend in NIR excess may 
be an indication of the youth of YSOs located in regions I and III as compared to those 
located near the center of Stock 8. To verify further the age sequence, we plot V/{V — I) 
and K/{H — K) CMDs for the three sub-regions in Fig. [221 which also suggests that the 
YSOs in regions I and III are relatively younger than those located near the cluster center. 

As mentioned earlier, several OB stars are located around the cluster Stock 8. We suspect 
that these sources belong to the first generation of stars formed in the IC 417 region, although 
the evidence is rather weak because of their scattered distribution. The age of the early- 
type stars in Stock 8 is estimated to be ^ 2 Myr, whereas the YSOs in the Stock 8 region 
indicate an average age of ~ 3 Myr (see Figs. fTOl and 1291) . However these two age estimates 
probe different physical mechanisms, viz., the former being a nuclear age and the latter a 
contraction age. So, we conclude that they are practically in agreement. We suspect that the 
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Figure 29. V/ {V — I) and K/ (H — K) CMDs for the Ha (triangles) and NIR excess sources (circles) lying in three sub-regions 
as shown in Fig. 1221 The isochrone of 2 Myr (continuous curve) by Girardi et al. (2004) and PMS isochrones (dashed curves) 
for ages 0.5, 5 and 10 Myr by Siess et al. (2000) corrected for cluster distance and reddening are also showrn. 



formation of stars in Stock 8 region was triggered by the first-generation stars mentioned 
above. The UV radiation from the OB stars, in particular the 08/09 stars (star numbers 
2 and 3 of Table 3), seem to have swept the pre-existing molecular cloud towards east of 
these stars. It formed a compressed gas layer which became unstable against self-gravity and 
collapsed to form second generation stars in the form of Stock 8 and surrounding YSOs in 
the region II. The shock/ionization front seems to be moving further to the east into the 
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low density gas and accumulating a denser layer of molecular gas, which we now observe as 
a hole in the MSX A-band emission region. 

The younger age of the YSOs scattered in region III is very strange. We suspect that they 
recently formed in scattered remnant clouds as seen in the Orion region (Ogura & Sugitani 
1998). A hint of the presence of such remnant clouds is noticed in the inset of Fig. [221 In 
this sense they may well be considered as stars of the third generation in the Stock 8 region. 

The narrow strip along the Nebulous Stream appears to be another active region (region 
I) of recent star formation. Fig. [22] shows that a few NIR-excess sources (11 stars) and 
Ha stars (13 stars) are distributed along the Nebulous Stream which also contains the CC 
14 cluster. The peaks of ^^CO and MSX A-band rioo contours along the Nebulous Stream 
almost coincide with the location of the embedded cluster. As we have discussed above, 
these stars show an average age of ~ 1 Myr. The question is whether the star formation here 
has been triggered by the UV radiation from the 0/B star in Stock 8 or from the O stars 
located further west. This is unlikely, because, as we have discussed in Sec. 10, the action 
of the UV radiation from these stars does not seem to reach this region. The morphology of 
the Nebulous Stream as well as the slightly-displaced distribution of the YSOs and CC 14 
to the south with respect to the Stream suggest that the pre-existing molecular gas might 
have been compressed by the shock front associated with the ionization front caused by the 
08 star located northwards of the Stream (cf. Section 11). So, the star formation activity 
along the Nebulous Stream is probably independent of that in the Stock 8 region. 



13 SUMMARY 

On the basis of a comprehensive multi-wavelength study of the IC 417 (Sh2-234) region we 
have made an attempt to determine the basic properties of the cluster Stock 8 as well as 
to study the star formation scenario in the region. Deep optical UBVIc and narrow-band 
[S II] , Ha photometric data, slitless spectroscopy along with archival data from the surveys 
such as 2MASS, MSX, IRAS and NVSS are used to understand the star formation in and 
around the cluster Stock 8. 

Reddening {E{B — V)) in the direction of the cluster is found to be varying between 
0.40 to 0.60 mag. The post-main-sequence age and distance of the cluster are found to be 
^ 2 Myr and 2.05 ±0.10 kpc respectively. Using 2MASS NIR two-colour diagram and grism 
survey for Ha emission stars, we identified candidate YSOs. A significant number of YSOs 
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(22 stars) are aligned along a Nebulous Stream eastwards of the cluster and an embedded 
cluster is located along the Stream at a distance of ~13 arcmin (7.8 pc) from the center of 
Stock 8. The mass of the YSOs lie in the range of ~ 0.1 - 3.0 Mq. The position of YSOs on 
the CMDs indicates that the majority of these stars in cluster Stock 8 have ages ~ 1 to 5 
Myr indicating star formation in the cluster may be non-coeval. 

The embedded cluster in the Nebulous Stream is found to be physically connected to 
the H II region. The YSOs located in the Nebulous Stream, embedded cluster and in the 
western region of Stock 8 have larger [H — K) excess in comparison to those located in the 
central region of Stock 8. The {H - K) excess, V/{V - I) and K/{H - K) CMDs indicate 
that these YSOs are younger than those located in the central region of Stock 8. The radio 
continuum, MSX, IRAS maps and the ratio of [S II] /Ha intensities indicate that the eastern 
region of the cluster is ionization bounded, whereas, the western region of the Stock 8 is 
density bounded. The morphology of radio emission and MSX A-band emission suggests 
that ionization/ shock front caused by the central ionization source of Stock 8 and by O 
type stars located in the western region of Stock 8, has not reached the Nebulous Stream. An 
08 star located ~ 9 pc away towards the north of the Nebulous Stream may be a probable 
source of the bright-rimmed Nebulous Stream. 

The slope of the mass function F in the mass range ~ 1.0 ^ M/Mq < 13.4 can be 
represented by —1.38 ± 0.12, which agrees well with Salpeter value (-1.35). In the mass 
range 0.3 ^ M/Mq < 1.0 the mass function is found to be shallower with F = —0.58 ± 0.23 
indicating a break in the slope of the MP at ~ IMq. The slope of the /^-band luminosity 
function for the cluster is found to be 0.31 ± 0.02 which is smaller than the average value 
(~ 0.4) obtained for young star clusters (Lada et al. 1991; Lada & Lada 1995; Lada & Lada 
2003). 
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